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Using lattice QCD we study the spectrum of low-lying fermion eigenmodes. According to the 
Banks-Casher relation, accumulation of the low-mode is responsible for the spontaneous break- 
ing of chiral symmetry in the QCD vacuum. On the lattice we use the overlap fermion formulation 
that preserves exact chiral symmetry. This is essential for the study of low-lying eigenmode dis- 
tributions. Through a detailed comparison with the expectations from chiral perturbation theory 
beyond the leading order, we confirm the senario of the spontaneous symmetry breaking and de- 
termine some of the low energy constants. We also discuss on other related physical quantities, 
which can be studied on the lattice with exact chiral symmetry. 
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1. Introduction 

Spontaneous breaking of chiral symmetry is the most fundamental property of the vacuum 
of Quantum Chromodynamics (QCD). Once we assume that the chiral symmetry is spontaneously 
broken, we can derive many important relations in the phenomenology of strong interaction, such as 
the GMOR relation, Goldberger-Treiman relation, and other soft pion theorems, which are written 
in the language of chiral effective theory. The problem of how and why the spontaneous chiral 
symmetry breaking occurs remains a difficult question due to the non-perturbative dynamics of 
QCD. 

Chiral symmetry of course plays a key role in the understanding of chiral symmetry breaking. 
In the flavor non-singlet sector of chiral symmetry, pion arises as the Nambu-Goldstone boson 
associated with the spontaneous symmetry breaking, while in the flavor-singlet sector the chiral 
symmetry is violated by the axial anomaly and is related to the topology of non-Abelian gauge 
theory. There are near-zero modes of quarks associated with the topological excitations; their 
accumulation in the vacuum leads to the symmetry breaking in the flavor non-singlet sector as 
indicated by the Banks-Casher relation. Therefore, the initial setup to study the chiral symmetry 
breaking should preserve both the flavor singlet and non-singlet chiral symmetries. 

Lattice QCD is the most promising approach to solve the low-energy dynamics of QCD, but 
there is a problem in realizing the chiral symmetry on the lattice. The conventional Wilson-type 
fermions violate the chiral symmetry at the action level, and the discrimination between the phys- 
ical effect of symmetry breaking and the lattice artifact cannot be done in a clear manner. On the 
other hand, the staggered fermions have a chiral symmetry but break the flavor symmetry. With 
these lattice fermions, the continuum limit has to be taken before analyzing the data with the con- 
tinuum chiral effective theory. 

Among other physical quantities, we are interested in extracting the chiral condensate {qq), 
which is an order parameter of the chiral symmetry breaking. This is not easy because the scalar 
density operator qq has a power divergence of the form niq/a^ as the cutoff \/a goes to infinity, 
hence the massless limit has to be taken to obtain physical result. (When the chiral symmetry is 
violated from the outset as in the Wilson-type fermions, the divergence is even stronger \/a^.) 
The condensate however vanishes in the massless limit, as far as the space-time volume is kept 
finite. Therefore, the proper order of the limits is to take the infinite volume limit first and then the 
massless limit, which is called the thermodynamical limit. Thus, the study of symmetry breaking 
is numerically so demanding, and some theoretical guidance is required to control the hmits. In 
QCD, the chiral perturbation theory (ChPT) provides such a theoretical framework. 

This work is an attempt to simulate the QCD vacuum on the lattice with exact chiral symmetry. 
We use the overlap fermion formulation [[I], ^, that exactly preserves chiral and flavor symmetries at 
finite lattice spacing and correctly reproduces the axial anomaly. We calculate the chiral condensate 
in various ways, which provides a good test of the chiral effective theory. In particular, we study 
the spectral density of the Dirac operator and compare it with ChPT at the next-to-leading order 
(NLO) of the chiral expansion. We also discuss a few other consequences of spontaneous symmetry 
breaking and related lattice calculations, which are made possible with exact chiral symmetry on 
the lattice. 
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Figure 1: Theoretical expectation from [Q] on the spectral function at L ^ 2 fm. The left and right panels 
correspond to the pion mass around 300 MeV and 100 MeV, respectively. The curves correspond to the full 
NLO calculation (solid curve, red), the leading order result in the e-expansion (dashed curve, blue), and the 
result in the infinite volume (dotted curve, black) at NLO of the conventional /7-expansion. 



These works have been done by the JLQCD and TWQCD collaborations. An overview of 
their recent physics results is found in [^]. 

2. Dirac spectrum and chiral symmetry breaking 

The chiral symmetry breaking is induced by an accumulation of low-lying eigenstates of 
quark-antiquark pair, as indicated by the Banks-Casher relation [||] 

I 

hm lim p(A =0) = -, (2.1) 

where p(A) denotes the eigenvalue density of the Dirac operator, p(A) = — Xk))- 

The expectation value (•••) represents an ensemble average and k labels the eigenvalues of the 
Dirac operator on a given gauge field background. On the right hand side of ( |2.1[ ), £ is the chiral 
condensate, £ = —{qq), evaluated in the massless quark limit. In the free theory, we expect a 
scaling p(A) ~ A^, which vanishes at A = 0, for a dimensional reason. The relation ( |2.1| ) implies 
that the spontaneous chiral symmetry breaking characterized by non-zero £ is related to the number 
of near-zero modes in a given volume after taking the thermodynamical limit. 

Based on ChPT, more detailed forms of p (A) at finite A, V and m have been obtained. This is 
achieved by evaluating the chiral condensate at imaginary value iX of valence quark mass, relying 
on the analytic continuation. The spectral function is thus obtained without explicitly treating the 
quark eigenstates. 

In this work, we use the most recent calculation by Damgaard and Fukaya [^], which is based 
on the conventional /^-expansion of ChPT but includes an integral over zero-momentum modes of 
pion field, hence it is also consistent with the £-expansion. At NLO, they provide a formula for 
p(A) at finite V and (non-degenerate) m. A typical example is shown in Figure |l| (left panel). The 
plot shows the spectral function p (A ) at a given volume V = x (3L) with L ~ 2 fm and pion 
mass niji ~ 300 MeV. Within the leading order (LO) of the e-expansion, the spectral function is 
given by the dashed curve, which is equivalent to what one obtains from the chiral random matrix 
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theory. It receives a finite volume correction at NLO and becomes the soUd curve. Starting from 
the second peak of the oscillating curve, the spectral function is suppressed due to the pion-loop 
effects. In the infinite volume limit, we expect a smoother dotted curve, which contains a chiral 
logarithm. 

The right panel of Figure |l| shows the spectral function in the £-regime. Here the lowest-lying 
eigenvalue is strongly suppressed by the fermion determinant as the quark mass is close to zero. 
The difference between LO and NLO in the £-expansion is less significant, since the system is in 
the £-regime. 

Once we could calculate the spectral density on a finite volume lattice, we can extract £ and 
1 /F^. Essentially, £ determines the height of the distribution and \/F^ the size of the NLO effects. 
The finite volume scaling should also be tested with more than one lattice volumes. 

3. Lattice analysis of the spectral density 

We use the lattice data obtained in the course of dynamical overlap fermion simulations by the 
JLQCD and TWQCD collaborations. The project mainly aimed at controlling the chiral extrapola- 
tion of physical quantities by realizing exact chiral and flavor symmetries in the lattice simulations. 
The use of the continuum ChPT is justified even at finite lattice spacings in contrast to other lattice 
fermion formulations, for which some modifications of ChPT with additional parameters is re- 
quired. With the exact chiral symmetry, there have been several new physics applications proposed 
from the project [^]. 

We use the overlap-Dirac operator D{Q) = [p /a)[\ +X /Vx^] with the Wilson kernel X = 
aDw — p [111 @]- This fermion formulation exactly preserves a modified version of the chiral sym- 
metry at finite lattice spacings [^]. The axial Ward-Takahashi identities are essentially the same as 
those in the continuum theory, and the index theorem is satisfied. 

The large-scale Monte Carlo simulations had been made feasible by restricting the Markov 
chain in a given topological charge, since the overlap-Dirac operator has a discontinuity on the 
boarder of the topological charge, and the numerical costs grows as to treat the change of 
topology. This is achieved by introducing unphysical heavy Wilson fermions in the simulation 
[^]. Fixing topology induces a finite volume effect for any physical quantities, but corrections 
are possible using the general formulae developed in |^ |T^. For the calculation of the spectral 
function considered in this work, the fixed topology is an advantage rather than a disadvantage, as 
the ChPT formulae are given for fixed topological sectors. 

At an early stage of the project, we performed simulations of two-flavor QCD at a lattice 
spacing a ~ 0.12 fm on a 16^ x 32 lattice. The simulation details are described in We took 
six values of sea quark mass in the range m.v/6 ~ nis, with rris the physical strange quark mass, to 
investigate the chiral extrapolation as discussed below. We also carried out a run in the £-regime 
with the quark mass around 3 MeV, in order to study the low-lying Dirac eigenvalue spectra, as 
discussed in the next section. 

We then extended the work to 2-1-1 -flavor simulations at a lattice spacing a ~ 0.11 fm on 
16^ X 48 and 24^ x 48 lattices. We cover a similar range of (degenerate) up and down quark masses 
as in the two-flavor runs while taking two values of strange quark mass near its physical value. A 
dedicated run in the £-regime has also been performed with the up and down quark masses around 
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Figure 2: Lattice results for the spectral function (histogram) together with the LO (dashed, blue) and NLO 
(soUd, red) curves. The left and right panels correspond to the pion mass around 300 MeV (/^-regime) and 
100 MeV (e-regime), respectively. 



3 MeV while keeping the strange quark mass near the physical value. A preliminary report of these 



runs is found in [ 1 1 ] 



At an early study of the eigenvalue spectrum, we use the distribution of the lowest-lying Dirac 
eigenvalue to extract £ in two-flavor QCD, by matching the distribution to the expectation from 
the chiral random matrix theory [ jl2| , |T^]. In this method the relation to ChPT is established only 
at LO of the £-expansion, so that the result may contain significant finite volume effect, which is 
the NLO effect. This was signalled already by a slight inconsistency between £ obtained from the 
lowest and from the second lowest eigenvalues. 

With the new formula ||5|], we can now consistently incorporate the NLO effects in the analysis. 
We use the 2-1-1 -flavor data of spectral function in a wider region of the eigenvalue A. Since the 
formula is valid also in the p-regime, we are able to include the /j-regime lattices into the analysis. 

Figure |^ shows the results for the spectral function in 2-1-1 -flavor QCD obtained on 16^ x 48 
lattices. The global topological charge Q is fixed to 0. The lattice data (histogram) are overlaid 
on the expectation from the NLO formula for both the /7-regime (left panel) and £-regime (right 
panel). We observe that the NLO formula nicely reproduces the shape of the lattice data. The 
curves are drawn by fixing the parameters (£ and F) with the data for an integrated spectrum at 
two representative values of A. Typically, one is taken near the top of the first peak, and the other 
is taken at A ~ 0.04, where the NLO effect is significant. 

We have further tested the NLO formula by extending the calculation to a larger lattice. The 
finite volume scaling is tested with the data on a 24^ x 48 lattice in the /j-regime. We find that the 
result for £ from this lattice is consistent with that on a smaller lattice after correcting the finite 
volume effect for both lattices. This implies that the lattice result scales as expected towards the 
thermodynamical limit. 

From this analysis, we obtain £ at six values of up and down quark masses m„j while keeping 
the strange quark mass at its physical value. Thus, we obtain an "effective" value of £ at each 
sea quark mass, L{mud,ms). We extrapolate the results to the chiral limit of mud using the NLO 



5 



Spontaneous chiral symmetry breaking on the lattice 



Shoji Hashimoto 



0.006 



0.005 - 




Figure 3: Chiral extrapolation of E(m„j,m5) to the limit of vanishing up and down quark masses. 



chiral expansion in two-flavor QCD [14] 



i;(m„d,m,) = i;(0,m,) 



1 



3Ml 



(3.1) 



as shown in Figure ^. (In the actual analysis, we use the formula including finite volume correc- 
tions.) 

Since we have the data in the £-regime, which is very close to the chiral limit, the chiral 
extrapolation is stable against the change of the fitting region. We also attempt to use the 2-1-1- 
flavor formula at NLO, and find the result in the chiral limit is consistent. 

After renormalizing to the MS scheme at 2 GeV using a non-perturbatively calculated renor- 
malization factor [|l|], we obtain a preliminary result I^^(0,m,;2 GeV) = [243(4)(+\^) MeV]^ 
where the second error represents an uncertainty due to the lattice scale. We also obtained F and 



Lg that are consistent with phenomenological analysis. For more details, see [16] 



4. Other consequences of spontaneous symmetry breaking 

The accumulation of the near-zero modes as seen in the spectral function is a direct mea- 
sure of the spontaneous chiral symmetry breaking, as its infinite volume limit corresponds to the 
Banks-Casher relation. On the other hand, there are many other physical quantities that reflect the 
spontaneous breaking of chiral symmetry. In the following, we describe a few of them calculated 
on the same set of lattice simulations. 

4.1 Topological susceptibility 

Topological susceptibility Xt = (2^)/^ characterizes how much topological excitations are 
occurring in the QCD vacuum. Although the definition involves the global topological charge Q, 
the susceptibility itself is a local quantity and could be determined even when Q is kept fixed [p^]. 
Namely, Xt can be extracted from topological charge density correlator as 

lim (mP(x)mP(0))e = f;^, - ^ + ••.')+ 0(^'-'"'-'W), (4.1) 
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Figure 4: Topological susceptibility Xt in 2 (circles) and 2+1 -flavor (triangles) QCD as a function of up and 
down quark masses. 



where mP{x) is a flavor-singlet pseudoscalar density operator related to the topological charge 
density through the axial-anomaly relation. When 2 = 0, this correlator approaches a negative 
constant at large separation \x\ after the excitation of the X]' meson saturates. This is intuitively 
understood as follows. Since the global topological charge is fixed to zero, if we find a positive 
topological charge excitation at a space-time point 0, then we have more chance to find a negative 
excitation at other points x to keep the total to be zero. 

This constant correlation is indeed observed in our work where we calculate the flavor-singlet 
pseudo-scalar correlator by maximally using the low-lying eigenmodes exactly calculated and 
stored on disks [17, [18]. We thus extract Xt at each sea quark mass. 

Sea quark mass dependence of Xt is plotted in Figure ^. We observe a good agreement with 
the expectation from the chiral effective theory at the leading order, = r/(l/m„ + l/nij + 



[]19[]. (m, is sent to infinity for the two-flavor case.) It provides another method to extract Z. Our 
result is [247(3) (2) MeV]^ (up to the error due to the scale setting), which is in agreement with the 
determination from the spectral function. 



4.2 Vacuum polarization functions 

The vacuum polarization function YI^j^'^\q^) defined through 

d\e''^''{0\J^{x)jl{Om = (g^q'-q^qM}\Q^)-q^qvnf\Q^) 



(4.2) 



for vector (7 = V) or axial- vector (7 = A) current is another probe of the spontaneous symmetry 
breaking. For example, the difference between vector and axial-vector channels is related to /;;: and 
Lio (one of the NLO low energy constants in ChPT) as 



fl = - lim 



-10 



hm TT— J 



Q 



n 



(1+0) 



n 



(1+0) 



(4.3) 
(4.4) 
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Figure 5: Difference of the vacuum polarization functions, Ily\Q^) — II^^\q^), multiplied by Q^. Data are 
shown for different quark masses; the chiral limit is also plotted by solid (linear extrapolation) and dashed 
(chiral logs) curves. 



which are called the Weinberg sum rules []20|]. Another interesting quantity is the electromagnetic 
mass difference of pion, which is expressed as 



Am, 



3 a, 



■EM 



4Kfi Jo 



n 



(1+0) 



(4.5) 



in the limit of massless pion. Since the combination {VV —AA) vanishes unless the chiral symmetry 
is broken, these quantities signal the chiral symmetry breaking. 

For the lattice calculation of these quantities, exact chiral symmetry is essential, since we have 
to extract a tiny difference between the vector and axial-vector channels. We extracted the dif- 



ference n''y\Q^)-n)^'{Q^) successfully, using the overlap fermion on the lattice 023], as shown 



r(l) 



in Figure |. The intercept and slope of Q'-[n'{^\Q^) -U^' {Q^)] at = correspond to and 
Lio, respectively, and the integral over gives Am^. The results clearly show that the sponta- 
neous symmetry breaking induces these physical quantities as expected. By further improving the 
numerical data especially in the low region, we will be able to precisely extract these quantities. 

Another use of the vacuum polarization function is an extraction of the strong coupling con- 
stant by matching the lattice data with the OPE expression in the perturbative regime ||23|]. 

4.3 Pion mass and decay constant 

ChPT is organized as an expansion in terms of small and p^, but the region of convergence 
of this chiral expansion is not known a priori. Using lattice QCD, one can test the expansion 
and identify the region of convergence. With the exact chiral symmetry, the test is conceptually 
cleanest, since no additional terms to describe the violation of chiral symmetry has to be introduced. 
(With other fermion formulations, this is not the case. The unknown correction terms are often 
simply ignored.) 

For the pion mass nijc and decay constant /;;:, the expansion is given as 



r(l), 



rrin 



IB 



1 „ 
1 + -xlnx + C3X + 0{x ) 



fn = / [l -Jclnx + C4X + 6>(x^)] , 



(4.6) 
(4.7) 
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where m-n and denote the quantities after the corrections while m and / are them at the lead- 
ing order. The expansions and ( p7| ) may be written in terms of either x = 2m^/{4nf)^, 
X = 2ml/{4nff, or <^ = 2ml/{4nf;c f (we use a notation of f„ = 131 MeV). They all give an 
equivalent description at this order, while the convergence behavior may depend on the expansion 
parameter. 

Figure ^ shows the comparison of different expansion parameters [24] in two-flavor QCD. 



The fit curves are obtained by fitting three lightest data points with the three expansion parameters, 
which provide equally precise description of the data in the region of the fit. If we look at the 
heavier quark mass region, however, it is clear that only the -expansion gives a reasonable function 
and others miss the data points largely. This clearly demonstrates that at least for these quantities 
the convergence of the chiral expansion is much better with the -parameter than with the other 
conventional choices. This is understood as an effect of resummation of the chiral expansion by 
the use of the "renormalized" quantities and /„. In fact, only with the -expansion we could fit 



the data including the kaon mass region with the next-to-next-to-leading (NNLO) formulae []24[]. 

Whether or not ChPT can be used for kaon is an important question and potentially has a 
strong impact on the light hadron phenomenology. To investigate this question, we are currently 



extending the analysis to 2h-1 -flavor QCD [g5|, |2 
5. Conclusions 

In this talk, I demonstrate that the scenario of the spontaneous chiral symmetry breaking in 
the QCD vacuum is confirmed by lattice QCD simulations with exact chiral symmetry. In the 
calculation of the Dirac operator eigenvalue spectrum, the expectations from ChPT is precisely 
tested including the NLO effects. 

With the use of exact chiral symmetry, new possibilities to extract physics from lattice QCD 
calculations have been opened. They include the precise calculation of the topological suscepti- 
bility, the vacuum polarization functions, and a theoretically clean test of chiral expansion in pion 
mass and decay constant. From the project, there are several other interesting calculations of var- 



ious physical quantities, such as Bk [27], meson correlators in the £-regime [M], nucleon a-term 



and strange quark content []29|], and pion form factors [ 
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Some of the applications discussed in this talk have been left without so much progress since 
the early days of QCD (or even before). At last, numerical simulation of lattice QCD has caught 
up theoretical conjectures made in 1960s and 70s, but now starting from the first-principles. 

I thank the members of JLQCD/TWQCD for fruitful collaborations. The author is supported 
in part by Grant-in-aid for Scientific Research (No. 21674002). 
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